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Abstract 

Mathematical models with time dependent parameters are of great interest in financial Mathematics 
because they capture real life scenarios in the financial market. In this study, via the Lie group technique, 
we analyse evolution-type equations with time dependent parameters and give the general symmetry 
structure of these equations. In addition, we illustrate this method by looking at an example of exotic 
options called the power options. Our model parameters are time dependent and the option gives a 
continuous yield dividend at different time intervals. We present new solutions, satisfying the boundary 
conditions, to this important problem. 
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1 Introduction 

Mathematical models of financial derivatives were arguably pioneered by Black and Scholes (B-S) pQ] and 
this led to the scientific legitimization of the options market. Many empirical tests have since shown that 
the price from this model is fairly close to the observed price [2J. A large number of financial derivatives are 
now modeled after the B-S equation with general form 

V t + XS 2 V SS + YSV s -ZV = 0, (1) 

where S is the underlying asset price, V(S, t) is the option price and X, Y, Z include important financial 
parameters like the stock volatility and interest rate. In reality though, these parameters change over time 
throughout the duration of the option. This necessitated the need to study equations that capture a scenario 
wherein the parameters are time dependent, viz. 

V t + X(t)S 2 V ss + Y(t)SV s - Z(t)V = 0. (2) 
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In our study we give the general symmetry structure of these financial derivatives (particularly when the 
model parameters are time dependent) and the route to their solution satisfying accompanying boundary 
conditions. We illustrate this method by looking at an exotic option called the power options. 

Exotic options are financial derivatives that have more complex features than commonly traded products. 
They are used to capture the increasing complexities of the financial market - especially the nonlinear 
dependence of the option payoff on the asset price - and are mostly traded over the counter to suit client 
needs. Power options are generally characterized by the nonlinear payoff structure of the underlying asset 
price raised to a certain power (f3) at expiration time. Previous studies by Cox and Rubenstein [3] looked 
at the case of squared power options and Tompkins [ 3 ] subsequently produced the complete solution to the 
problem of the power parameter as any given integer value. Esser [5] produced solutions to the problem 
when the asset price is raised to any real valued number. Okelola et al [6| then looked at the realistic case 
of the model parameters being time dependent. 

We note that the PDE ([2]) is linear and a study by Johnpillai and Mahomed [8] has found criteria for 
parabolic linear homogeneous equations to be equivalent to the classical heat equation. They showed that 
if the relative invariant of the equation is zero, it can then be concluded that the equation is equivalent to 
the heat equation. It should however be noted that for the problem we look at in this paper, it will not 
be a straightforward matter to transform the solutions and determine the behavior of equation ([2]) from 
that of the heat equation. In addition, the solution to equation <[2|) is required to satisfy certain boundary 
conditions and this demand cannot be forced on the heat equation in a direct manner. Due to these reasons, 
we decided to examine the PDE apriori. 

2 Lie group approach 

The Lie group technique utilizes the notion of symmetry to obtain solutions to DEs in an algorithmic 
manner. Consider an nth order system of DEs 

~d(t,y 7 dy,...,d n y) = 0, d = l,...,m (3) 

where t = t 1 ,... ,t e are the independent variables, y = y 1 ,... ,yf are the dependent variables and d l is the 
l- th partial derivative of y with respect to t (for 1 = 1,... ,n). An infinitesimal generator 

X = C(t, y)d t i + rf(t,y)dyp, (4) 

of the one-parameter Lie group of transformations 

(t*) 1 = f l (t,y,e), (5a) 
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(y*r = g^t,y,e), 


(5b) 


is said to be point symmetry of equation Q if and only if 

X (n) H d (f,y,dy,... ,d n y) = 0, (6) 

where X ^ is the n th prolongation of X [7]. 

A local group of transformations Z (acting on M C T x Y = M e x WJ ) is said to be a symmetry group of 
S if every element z € Z transforms solutions of S to other solutions of S. On reduction of S via Z, every 
solution of the new system ( E/Z)d(a , b , db ,... , d n b ) gives rise to a Z-invariant solution to S. This technique 
has been used to great effect in Financial Mathematics mm- 


3 Symmetry structure of the problem 


The PDE © admits the six parameter symmetry [11] 


G — 9(t)dt + I 7 (t) + 


log (S)9(t) , log (S)X(t)d(t)\ 


+ 


Sd 8+ \V{S,t)a(t) 


+ f Z(t) log (S)9(t) + Z(t) log(S)0(t) - log(5)d(t) + [Y(t) log 2 (S') - X(t) 2 log 2 (S) 
— 2X(t) log(S)j ki \ V(S,t) / \Y{t)-X{t)) W, 


( 7 ) 


subject to the functions 9(t),a(t ) and 7 (t) satisfying 


X(t) 2 9{t) - 8X(t)h - X{t) 2 9[t) + X(t)X{t)9(t) + X{t)X(t)9(t) = 0, (8) 

- 3X(t) 2 Y{t)9(t) - 2 X{t) 2 Y{t)9{t) + 2X(tf^(t) + 3 X(t)X(t)Y{t)9{t) + 2 X(t)X(t)Y(t)9(t) 

+ 3X(t)X(t)Y(t)9(t) - 2X(t)X(t)j(t) - 3X(t) 2 Y(t)9t = 0, (9) 

2X(t) 2 Y(t)9(t) - AX(t) 2 Z(t)9{t) - X(t) 2 X(t)9(t) + 2X(t) 2 Y{t)9(t) - 4X{t) 2 Z{t)9(t) - 2X{t) 2 ^ 

+ 4 X{tfa{t) - X{tf9(t ) + 8X(t) 3 h + X(t)Y(t) 2 6{t) - X{t)Y(t) 2 0{t) - 2X(t)Y(t)Y(t)6(t) 

+ 2X(t)Y(t)i(t) = 0. (10) 


Since the equation is linear, it also admits an infinite-dimensional ‘solution’ symmetry |7i]. The symmetry 
structure in ([7]) and the accompanying equations in the system ([8]) (flOll is applicable to all evolution-type 
equations of the form ([2]) . We will show its applicability by considering an example of exotic options. 
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3.1 Solution to the power option problem 


We now consider the power option that pays out a continuous yield dividend, y(t), at different time intervals. 
To derive the PDE modeling this option, we assume that the stock pays a dividend at a continuous rate 
proportional to the value of the stock. The dividend paid per unit time will thus be y(t)S and the dividend 
paid in a short time interval will be y(t)Sdt. The pricing kernel is characterized by the Brownian motion 
model 


d'i p 

0 


—r(t)dt 


P0) ~ r (t) 

a(t) 


d W(t), 


( 11 ) 


and the stochastic process for S 13 is given by 


dS& = [(//(f) - y{t))/3 + |/3(/3 - 1)] S?dt + a{t)pS^dW(t), 


( 12 ) 


where r(f) is the stock interest rate, a{t) its volatility, //(f) the instantaneous expected return and W(t) is 
the Weiner process. By Ito’s lemma, the price of the power derivative will be 

dV = [^a{tffi 2 S^V sPsP + ^ + [(//(f) - y(t))P + - l)a(f) 2 ] S^V s ^j dt + a(t)f3S^V s ,W(t). (13) 

On solving for if: and V in equations (111|) and (1131) respectively, their product (since it is a martingale the 
drift term must be zero) gives the PDE 
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V t + ^( t) 2 p 2 S 2 V ss + (r(f) - y(t))/3 + ^0(0 - l)a(t) z SV s - r{t)V = 0, 


1 


(14) 


where the payoff is given by V ( S , T ) = max{T(S ,/3 — A'), 0} at the boundary t = T. The option has a payoff 
V(S,T ) = 0 when S 13 < K and a payoff V(S,T ) = S 13 — K when > K. The strike price is K and T 
determines whether the option is a call or put. 


On substituting the representations of X(t),Y(t) and Z(t) of equation (fill) into the system © - (fTOl) and 
solving, we obtain 


m = 


a(ty 
7 (t) = 72 + 


02 + j cr(t) 2 dt(j)i + 2k/3 2 j cr(t) 2 dt^ 

J dt 


<r(t) 7i + 0 


+ 


u(f) 2 ( 

3y(t)&(t) 4r(t)a(t) 2 4 y(t)a(t) 2 


3 r , 3m(f) 3y(t)a(t) 

-Y + ~Yt) -~ + 3r(i)a(t) 


a(t) 


a{ty 


+ 


a(tf 


+ y(t) + 


cr(f) 

2r(f)<r(f) 2y(t)d(t) 


a{t) 


u(f) 


m 


m 


df, 


a(t) = / { [ r(t)a(t) 3 - 2r(t)y(t) + r(t) 2 a{t) + y(t) 2 a(t) + y(t)(r(t) 3 + 


( tf\ m 

4 J 2a(t) 2 


(15) 


(16) 
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+ ( y(t) + ^- r(t)] - cr(t) 2 /3 2 Aq + ( 1 - r(t) 2 d(t) - y(tj 2 a(t ) 


2 ' K ~'J cr(f) 2 /3 ” v “' " \ 2 

+ 2 r(t)y(t)a(t) + r(t)r{t)a{t) - r(t)y(t)a(t ) + y(t)y(t)a(t) - y(t)r(t)a(t ) 
y(t)a(t) 3 cr(t)a(t) A 


+ 


+ 




| 0(f) J-di + cci, 


(17) 


2 ' 4 

where 0i, 02 , 71,72 and oq are constants of integration. 

3.1.1 Solution to the power option problem with payoff V(S,T) = 0 

Application of the general symmetry ([7J) to the terminal condition will demand that 0(T) = 0. For this to 
be true, the parameters 6 \, 02 and Aq must be zero. This requirement reduces the 6-paranreter symmetry to 
a 3-parameter symmetry. This new symmetry has invariants P and u resulting in 

( log(S) 


V = P{u) exp 


\ 7 (t)P 




where u = t. We use these to reduce the power option PDE to an ODE with solution 


P = Ci ( cosh ^uAi(t) j — sinh ^uAq(t) j ^(A)^ 


(18) 


(19) 


where the A,’s are known functions of time. The solution of the PDE satisfying the terminal condition will 
thus be 


V = C\ ^cosh (uAi(t) jA 2 (t) S j — sinh ^udi(t) j Aq(A)) 


i 

7(0/3 


s 


ai/3- 


71 log(S) 
2/3 


+71 / ( r(t)-y(t)-A Vdt 


( 20 ) 


where the constant C\ is determined by the requirement that V(S,T) = 0. 

3.1.2 Solution to the power option problem with payoff V(S,T) = S' 3 — K 

For the condition V(S,T) = S 13 — K at t = T, we follow the same steps as for the terminal condition. The 
general symmetry in equation (J7]) reduces to a one-parameter symmetry with the constants 0i, 7i, 72 , aq and 
k\ all equal to zero. The function 0(f) is redefined as 


e\t) = 


0 2 


1 — / cr(f) 2 df 


cr(f) 2 df 


( 21 ) 


t=T 


with obvious modifications for 7*(f) and a l {t ) (see (1161) and (1171) h The solution will then have the form 


V = P 1 (u) exp 


f ^(t) _ 

\0 1 (t) 2r(t)p-(5a(t) 2 -2/3y(t)[6 1 (t) 6\t) 


2 log (5) 


d 1 ^) r{t)6 l {t) 


-r(t) 


>dt, , 


( 22 ) 
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where 


u = 


/ 


T 1 ^) _ dt _ log(5) 


a(tw(tw *m(ty 

These reduce the power option PDE to a second order ODE with solution 


(23) 


P 1 = cosh ( 2 a(t) + ub(t ) — ‘ ^ ^ ) — sinh ( 2 a(t) + ub(t ) — ^ ^ 


b(t) 


b(t) 


ClH -i+ c M+ a(t)c{t) 


b{t ) 3 b(t) 


2 ’ 


y/2b(t) 2 


/ b(t ) 3 + a(t)b(t)c(t) — c(t) 2 1 

^ a(t)b[t ) + ub(t) 2 — 2 c(t) 2 ^ \ 

y 2b(tf ’2’ 

2b(ty JJ 


(24) 


where H and \F\ are the Hermite and Hypergeometric functions respectively. As before a(t),b(t) and c(t) 
are known functions of time. The solution to the problem with payoff V ( S , T) = S'' 3 — K will thus be 


V (5, T) = exp 




ClH - 1+ c M+ a{t)c(t) 


b(t ) 3 &(*) 


2 ’ 


i(t)b(t) + ub(t) 2 — 2c(t)\ (b(t) 3 + a(t)b(t)c(t) — c(t) 2 1 y a (t)b{t) + ub(t) 2c(t) 

/ 21 V ’ 2 ’ 26^p 


'I 


a(t ) — r(t)6(t) — f{t)6(t) 


-df 


5 J 0{2r(t)P - (3a{t) 2 - 2f3y(t)) 

The constants C\. C '2 are determined from the requirement that V(S,T ) = S' 3 — I\. 

4 Conclusion and remarks 


(25) 


We have shown that all evolution equations of the form (f2j) . regardless of the nature of the coefficients, will 
have the 6-paranreter symmetry structure (0 whose solution route will follow the same steps outlined in 
this paper. As a result, we can find group invariant solutions to a large class of financially relevant DEs. 

We introduced a new PDE modelling power options with time dependent parameters and paying out a 
continuous dividend at different time intervals. Our general result was then contextualized to this equation. 
New solutions satisfying the boundary conditions were provided for the first time. 

We believe that this approach will be useful in Financial Mathematics. 

Acknowledgements: We thank the University of KwaZulu-Natal and the National Research Foundation 
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